Macroscopic Quantum Dynamics of a Free Domain Wall in a 

Ferromagnet 

0\ • Junya Shibata^ and Shin Takagi^^ 

Department of Physics, Tohoku University, Sendai 980-8578, Japan 

(February 1, 2008) 



0^ 



o 

in 



o 



o 



1 



o 



% 



Abstract 



A. . We study macroscopic quantum dynamics of a free domain wall in a quasi- 

1^ , one-dimensional ferromagnet by use of the spin-coherent-state path integral 

in discrete-time formalism. Transition amplitudes between typical states are 
't^ I quantitatively discussed by use of stationary- action approximation with re- 

spect to collective degrees of freedom representing the center position and the 
^_^ chirality of the domain wall. It is shown that the chirality may be loosely said 

pj . to be canonically conjugate to the center position; the latter moves with a 

O \ speed depending on the former. It is clarified under what condition the center 

position can be regarded as an effective free-particle position, which exhibits 
the phenomenon of wave-packet spreading. We demonstrate, however, that 
^ ■ in some case the non-linear character of the spin leads to such a dramatic 

phenomenon of a non-spreading wave packet as to completely invalidate the 
Q \ free-particle analogy. In the course of the discussion, we also point out various 

^^ ' difficulties associated with the continuous-time formalism. 
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I. INTRODUCTION 

Recent developments in low-temperature measurement techniques and the so-called 
nanostructure technology enable us to study low- dimensional magnetism in mesoscopic mag- 
netic systems. Among others, dynamics of a domain wall in a ferromagnet has attracted 
much attention both theoretically and experimentally, because it is expected to exhibit 
quantum-mechanical aspects at sufficiently low temperatures.au A domain wall contains a 
(semi-)macroscopic number of spins, its width being typically 10 ~ 1000 A. Hence, if its 
quantum-mechanical behavior was found, it would be an evidence of macroscopic quantum 
phenomena (MQP). To list just a few of the theoretical works about possible MQP involv- 
ing such a domain wall: depinning of a domain wall via macroscooic quantum tunneling 
(MQT),a0 coherent tunneling though a periodic pinning potential,0 macroscopic quantum 
coherence (MQC) of the chirality.Erilj 

One of the standard procedures to discuss MQP begins by deriving an effective action 
in terms of those collective degrees of freedom which directly describe the tunneling in 
question. In the case of the magnetic domain wall, the relevant collective degrees of freedom 
are the center position and/or the chirality of the wall. Existing works in the literature 
then treat the effective action in the Caldeira-Leggett schemes to evaluate the tunneling 
rate. However, as emphasized by Leggett,li^ one should probe a quantum-mechanical time 
evolution to check whether MQP (especially, MQC) have occurred. Hence, what is needed 
on the theoretical side is to evaluate not only tunneling rates but also relevant transition 
amplitudes. 

As a technique to evaluate a transition amplitude, the spin-coherent-state path integral 
in continuous-time formalismllJ is frequently used. However, as noted by some worker s,I1§lJ 
this formalism has some fundamental difficulties even at the level of a single-spin system. 
The nature and implication of the difficulties have recently been examined in detail,li3 where 
it has been pointed out, among others, that the information on the initial and the final 
states fails to be retained in the transition amplitude in question, and that, when combined 
with the stationary- action approximation, the stationary value of the action is not always 
given correctly and the fluctuation integral diverges. Hence, at present, there are no reliable 
results for transition amplitudes. One of the interesting predictions made by the continuous- 
time formalism is an interference effect (the so-called spin-parity effect); the behavior of the 
domain wall is predicted to depend dramatically on whether the magnitude S of each of the 
individual spins is an integer or a half-integer. The origin of the effect has been ascribed to 
the so-called Berry-phase term appearing in the effective action.u However, the information 
of the initial and the final states are essential for interference effects. Hence, the purported 
interference effect need be re-considered. 

In this paper, we present a ffist step to clarify these problems by use of the discrete-time 
formalism of the spin-coherent-state path integral. We focus upon a free domain wall and 
evaluate real-time transition amplitudes. This is the first case of an unambiguous evaluation 
of such amplitudes via the spin-coherent-state path integral as applied to an interacting 
many-spin system. 

The paper is organized as follows. In Sec. II we present the model Hamiltonian to be 
treated, namely that consisting of the Heisenberg exchange and the anisotropy energies, and 
formulate transition amplitudes between spin coherent states by use of the spin-coherent- 



state path integral in the discrete-time formahsm. In Sec. Ill we introduce a domain wall 
together with the collective degrees of freedom for its center position and the chirality, re- 
spectively. This section is also devoted to the derivation of the effective action of the domain 
wall. We also point out some problems associated with the continuous-time treatment of 
the effective action. Sec. IV evaluates transition amplitudes for a free domain wall in the 
stationary-action approximation including the effects of fluctuations. In the course of the 
evaluation, we note the conjugate relation between the center position and the chirality. We 
also locate those terms which can induce interference effects. In Sec. V we compute transi- 
tion probabilities between typical states and compare the quantum dynamics of the domain 
wall with that of a free particle. This comparison allows us, among others, to identify the 
"effective mass" of the domain wall. In Sec. VI it is explicitly shown that the continuous- 
time formalism leads to a wrong transition probability for a free domain wall. We conclude 
with a speculation on the possibility of MQT and/or MQC involving a domain wall. 

II. MODEL 

We consider a ferromagnet consisting of a spin S of magnitude S at each site in a quasi- 
one-dimensional cubic crystal (a linear chain) of lattice constant a. The magnet is assumed 
to have an easy axis and a hard axis in the z and the x directions, respectively. Accordingly, 
we adopt the Hamiltonian 

H = -JY.S^- S, - - Y.{KSl - K^l), (2.1) 

where the index i or j represents a lattice point, < i,j > denotes a nearest-neighbor pair, 
A'^^: is the total number of lattice points, and J is the exchange coupling constant, and K 
and K_j_ are longitudinal and transverse anisotropy constants; J, K, and K_j_ are all positive. 
Since we are interested in those transition amplitudes which are appropriate to describe 
quantum mechanical motion of a domain wall, we introduce a spin- coherent statetB at each 
site, which is suited for a vector picture of spin. By use of the eigenstate \S) of Sz associated 
with the eigenvalue S, the spin coherent state is defined by 

|n):=(l + |er)-^exp(eS_)|5), (2.2) 

where n is a unit vector {rix = sin 6 cos 0, Uy = sin 6 sin 0, n^ = cos 6) with the complex 
number i^ being its Riemann projection: 

^ = e'HaJ-, r = e-^nan^. (2.3) 

These states form an overcomplete set and possess, among others, the following properties: 

{n\S\n) = Sn, (2.4a) 

(n'|(e-5)» 



{{n'\{e-Sy\n)): 



{n'\n) 
l-^)(((n1e-5|n)))2 + |, (2.4b) 



where e is an arbitrary unit vector. Hereafter we work with the ^-representation, and denote 
a state of the system as 

Nl 

\0^\^i,^2,---,U)--=<S)\Q^ (2-5) 

j 

where |^j)(= |^j)) is a spin coherent state at the site j. The transition amphtude between 
the initial state |^i) and the final state |^f) can be expressed as a spin- coherent- state path 
integral in the real discrete-time formalism by the standard procedure of the repeated use of 
the resolution of unity (see, e.g., Ref. 17 on which the present notation is based): 

(CFle-^^^/^lei) = ji^n^y n n^/^(0(^)>e;(^))exp [-^S[i\i\) , (2.6) 



n=l 



J 



where A^ = T/e, e is an infinitesimal time interval, n represents discrete time, and the 
integration measure is 

_ 2S+1 dQin)d^*in) dUn)di*{n) ^ d^^,{n)dQ^,{n) 

The action iS[^*,^] consists of two parts, 5'^[^*,^] and iS'^[.^*,.^], which are to be called the 
canonical term and the dynamical term, respectively. They take the following forms: 

S[C,^]=S^[C,^]+S^[C,^], (2.8a) 

„• N Nl 

T'5ir,e]:=EEMe,(n)|e,(n-l)) 

" n=l j 



n=l 



tS'[C, e] := -T E ^m*in),an - 1)), (2.8c) 

'^ '^ n=l 

H{c,v)--={m\v)), m^^h am^^F. (2.8d) 

Here, we emphasize that the integration variables are {C,*{n),C,{n)\n = 1,2, ..., A^ — 1}; ^(0) 
and ^(A^) are fixed complex numbers. In passing, note that 

Un)Un - 1)) = (cos^ ^ cos^ ^^^ + sin^ ^ sin^ M^e'^iMnyMn-i))^ , (2.9) 

which is a 27r-periodic function of the phase difference 4>j{n) — (f)j{n — 1). 

At this stage, if one regarded all the differences \C,j{n) —^j{n — 1)| as small in some sense, 
expanded the action iS[.^*,.^] and went over to the continuous-time formalism, one would 
obtain the following form which had been used in most of the literature including Refs. 8 
and 10: 



{^F\e-'^^/%) ~ / V^VC exp {^ {sue, ^] + 5,^ [C, ^l) } , (2.10a) 

'j-SLlC,^] := -l[dtH{C{t),m\ (2.10c) 

where V^VC,* represents a symbolic measure in the continuous-time formahsm, However, as 
pointed out in our previous paper ,0 this formahsm has various difficuhies (see also the next 
section). For this reason, we proceed to consider the transition amplitude in the discrete-time 
formalism ( p.8| ). 

We shah be interested in those spin configurations whose scale of spatial variation is 
much larger than the lattice constant a. Accordingly, we take the spatial continuum limit 
in (lObl) and (gl^: 



' S^\r fl - 5 V f'^" ^^ In (l + r(^.^)e(a:,n-l))^ 

^ ^^ '^J - \^,y-./2T'"(i + ie(x,n)p)(i + ie(x,n-i)p)' ^'-^'"^ 



n{C{x),r]{x)) :-- 



lS%\i] = -j-j:^eJ_^^^-n{C{x,n),ax,n-l)), (2.11b) 



:i+e{xMx)y 2^^^^-r(x)5..,(a:) 

-f {(5-^)(l-r(a:)r/(x))^ + i 

aK r/„ 1\ _. o 1 
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{s--){i-c{xHx)y + 
+^{{s~l){c{x)+v{x)r + 



(2.11c) 



where L is the length of the linear chain, J = J a?, and a = K^/K. In this paper we consider 
the case of a weak transverse anisotropy a ^ 1, and study the dynamics of a domain wall 
to the first order in a. 



III. EFFECTIVE ACTION FOR A DOMAIN WALL 

A. Kink configuration 

We begin by finding a domain wall configuration. It is determined by one of the static 
solutions {^^(a;),,^'^(a;)} of the action iS[^*,^]. They satisfy the following equations up to 
O(a0): 

A-ferM- ™gyf }-i-ffMeM^O. (3.1a) 

I l + r(a;)^'(a;) J l + ^"(x)^'^(x) 

1''^^ l + e(a;)e(a;) J l + e(2;)f(a;) 



where A^ = JS/K{S — 1/2). An obvious solution is the "vacuum" solution representing the 
uniform configuration in which the spins are either all parallel or all anti-parallel to the z 
direction. The other solution is the "kink" solution representing a domain- wall configuration 
in which the spins at x ~ -|-oo are parallel to the z direction, the spins at x ~ — oo are anti- 
parallel to the z direction, and there is a transition region (i.e., a domain wall) of width 
A; 

r (x) = exp (-^^ + #o) , f (x) = exp (-^^ - ^<Po) , (3.2) 

where Q and 0o are arbitrary real constants. Q is the center position of the domain wall, 
and 00 is a quantitative measure of the chirality of the domain wall with respect to the x 
axis (Fig. 1); the wall is maximally right-handed if 0o = 7r/2 and maximally left-handed if 
00 = ~7r/2, while it has no chirality if 0o = 0. The range of 0o is chosen as — tt < 0o < tt, 
with 00 = vr and 0o = — vr representing the same situation. {^'^(— a;),^'^(— x)} is also a 
solution representing a domain-wall configuration. However, this as well as the vacuum 
solution belongs to a sector different from that of (|3.2| ). Since a transition between different 
sectors are forbiddenjl^l it is sufficient to consider only the sector ( |3.2|) for the purpose of 
studying the dynamics of a domain wall. 



B. Collective degrees of freedom 

Study of the domain- wall dynamics is facilitated by introducing relevant collective degrees 
of freedom. We note two kinds of invariance possessed by ( p.ll) . One is the translation 
invariance in the x direction. The other is the rotation invariance around the z axis. These 
invariances are embodied by the arbitrariness in the choice of Q and cbn. respectively, in 
( p.2[ ). Hence we elevate them to dynamical variables Q{n) and 0o(?7.).&t3'EJ To deal with 
these two dynamical variables (collective degrees of freedom), it is convenient to define 

z{n) := q{n) + i(f)o(n), z*{n) := q(n) - i(f)o{n), (3.3) 

q{n) = Q{n)/X, n = 1,2, ..., N - 1. 

By use of these variables, original integration variables ^(x, n) and ^*(x, n) may be decom- 
posed into the domain- wall configuration and the environment around it: 

,^(x, n) = Ci^'i -2(^)) + v{^^ ^1 -2(^))) (3.4a) 

C{x,n)=C{x;z*{n)) + 7]*{x,n;z*{n)), (3.4b) 

where 

C{x;z{n)) ■=exp{-x/X + z{n)), Cix; z*{n)) := exp {-x/X + z*{n)) . (3.5) 

Hereafter we consider transition amplitudes between the following domain-wall states: 

1^/3) = \zp) := (g) If (ja; z^)), (3 = I,F, (3.6) 



where zi represents the center position qi and the chirahty 0i of the domain wall in the initial 
state, and zp those in the final state. At both ends of the discrete time (n = or n = A^), 
we define 



z{0) = Zi := qi + i(j)i, z{N) = Zf := gp + #f, 

ri{x, 0; z{0)) = ri*{x, N; z*{N)) = 0. 



Putting Eqs. 



into the action (|2.11|), we obtain up to 0{a) 



<S[Cj^] = S^[z*, z] + terms involving the environment 7], 

S'[z*,z] = S'^'lz* , z] + S'^'lz* , z], 



? / (JiiAj 111 



^5«|.-..]. ^ 



N 



n _|_ (i-x+z*(n)+z{n)\^-^ I ^—x+z*{n—l)+z{n—l)\ ' 



Is^'[z\z]:=-'-E^^J:^ 

n=l 



h^ 



a 



1 + - {1 + cosh(2*(n) - z{n - 1))} 



(3.7a) 
(3.7b) 



(3.8a) 
(3.8b) 

(3.8c) 
(3.8d) 



where A'dw = A/a is the number of spins in the domain wall, and -Edw = 2AdwA'5'(S' — 1/2) 
is the kink energy. The zero point of energy has been adjusted in (|3.8d|). In this paper. 



we consider only 5^[2;*,z] which is expected to make the most dominant contribution to 
transition amplitudes. The infiuence of the environment shall be discussed in a separate 
paper. 

Expression ( |3.8c|) can be reduced to a simpler form (see Appendix A for the details of 
the derivation by use of dilogarithnv^): 



h 



S''[z*,z]/Nny,S 



1 ^ 

-y 



(g(n) - qin - l)f - RiMn) - M^ - 1)) 
-i{2{q{n) + qin - 1)) + 2L/A}/(0o(r^) - 0o(^ - 1)) 



where both R{(f)) and J(0) are 27r-periodic functions (see. Fig. 2) such that 

R{4>) = 02 : 101 < vr, 

/(0) = : — TT < < TT. 



(3.9) 



(3.10a) 
(3.10b) 



This periodicity follows inevitably from the periodicity of ( |3.8c| ). In general this periodicity 
need be respected in performing the integration with respect to {0o(?t-)}. However, in special 
circumstances when the range of 0o(^) — (poin — 1) can be restricted to [— vr, vr] for all n, the 
above action may be rewritten as 



-S''[z*,z]/Nj,^S 



E - o (('^H - <li^ - 1))' + (Mn) - Mn - I)?) - tiqin)Mn - 1) - Mn)qin - 1)) 

n=l *- ^ 



-i {2(0F - 0i)-^/A + (gF0F - gi0i)} 
^ r 1 

= Y^ < — {z*{n)z{n) + z*{n - l)z{n - 1)) + z*{n)z{n - 1) 

n=l *- ^ 

-I {2(0F - (t^i)L/\ + (gp^F - gi0i)} . (3.11) 

This expression, except for the last constant term, formally coincides with the corresponding 
action appearing in the (boson-) coherent-state path integral with a single degree of freedom. 
The last term can be neglected because it is a constant phase, which does not affect any- 
physical quantity. 

C. Effective action for a domain wall in continuous-time formalism 

Let us comment on the continuous-time treatment of a domain wall and the associated 
problems. 

If one started from (|2.10|) , one would obtain the continuous-time counterpart of (|3.^ ) as 



Jo 



n := -^^(r(t)z(t) - z*it)zit)) - Ebw 



l + ^{l + cosh(z*(t)-^(t))} 



= 2^^ <5(t)0o(t) - 2A/'Dw^^ficos^0o(t), (3.12a) 

A 

0.|(s-i)a. (3.12b) 

In the last expression, constant terms including those appearing as a result of partial in- 
tegration have been neglected. If the transverse anisotropy is strong in the sense that 
A^DW'S'fi ^ 1) 00 (^) Kiay be restricted to a region near ±7r/2; 

0o(t) =C77r/2 + ¥;(t), C = ±l, |¥^(t)|<l. (3.13) 

Substituting this into ( |3.12a| ), one would find 

£con ^ CAQ{t) + '^AQ{t)Lp{t) - 2NnwfiSn^\t), (3.14) 

TT 

where A = N-D^^hSn/X. Provided that the path-integration measure is independent of ip, 
Gaussian integration with respect to (f would then lead to the following effective action for 
Q: 



'j-S:L[Q] = 11 dt [cAQit) + -fQ'it)] , (3.15) 

where Md is the Doring mass: 



Mt,= 



2h'^ 



a 



2 



KS 1 , , 

(3.16) 



^ J(^-l/2)i^±' 



Accordingly, one might expect the following expression for the transition amplitude: 

{z^\e-^^^/^\z,) ~ /l^gexp (^5l[g]) . (3.17) 

This corresponds to the result of Braun and Loss.S At this stage, it has been concluded that 
the center position of the domain wall behaves as a free particle with a possible modification 
due to interference effects induced by the first term CAQ{t) of the effective action, which is 
often called the " Berry-phase term" . 

If one is interested in the quantum depinning of the domain wall, a pinning potential 
is to be added to Lagrangian ( |3.14| ). If the effect of the transverse anisotropy is relatively 
larger than that of the pinning potential, the effective action ( p.l5|) is simply augmented by 



the pinning potential. MQT of the center position has been discussed on the basis of this 
action. On the other hand, if the pinning effect is the larger, one could carry out Gaussian 
integration with respect to Q{t) and obtain an effective action governing (/>o(i), with which 
MQC of the chirality has been discussed.llj 

However, the whole series of the above-quoted arguments, which are based on the form of 
the derived effective actions regardless of the path-integration measure, are at best heuristic 
and their validity is rather dubious. In the literature it has been tacitly assumed that the 
right-hand side of ( p.l7| ) is a Feynman kernel with respect to Q. (Recall the following point: 



In order to obtain a transition amplitude between physical (i.e., normalizable) states from 
a Feynman kernel, the latter has to be multiplied by the initial and the final wave functions 
and integrated over Q-p and Qi.) The starting point of the whole arguments, on the other 
hand, is the left-hand side of ( |3.17D , which is a transition amplitude between physical states. 
Its path-integral structure is different from that for a Feynman kernel. One should not be 
misled by the apparent form of the effective action (|3.15|) . Though integration over ip{t) 
may be carried out in principle, resulting action can not be a Feynman kernel. Indeed, 



as we shall illustrate in Sec. VI, the right-hand side of (|3.17|) as interpreted as a Feynman 
kernel does not give a correct transition amplitude. By the same token, interference effects 
predicted on the basis of the "Berry-phase term" CAQ{t) need be re-examined: as shown in 
the following section, since the continuous-time treatment neglects many other terms which 
can contribute to interference effects, there is no reason that only the "Berry-phase term" 
should be kept. 

IV. TRANSITION AMPLITUDE IN STATIONARY- ACTION APPROXIMATION 

In this section, we evaluate transition amplitudes by means of the stationary-action 
approximation, discuss the conjugate relation between the center position and the chirality 
of the domain wall, and locate those terms which can induce interference effects. 



A. Stationary-action path 

Let {z^ {n) , z^ {n)} be the stationary-action path, namely, the stationary point of the 
action S''[z*,z\: 



dS'[z*,z] 



dz*{n) 

dS^'lz^z] 



dz{n 



0, 
^0, 



n 



1,...,N-1, 



n = 2,...,N, 



(4.1a) 
(4.1b) 



where the symbol |sindicates the replacement {z*,z) -^ {z^,z^) after differentiation. It is 
convenient to defineO 



z'{0) ■=zi = qi + i^ 



z^'i^N) := Zp = qF-i^p. 



(4.2) 



Let us work with (|3.11|) instead of (|3.9| ). This procedure will be justified a posteriori. Then, 
the above set of equations take the form: 



z''{n) - z''{n - 1) = -ieVLsmh.{z''{n) - z''{n - 1)}, 
z''{n) — z''{n — 1) = —ieVt?,mh.{z''{n) — z''{n — 1)}, 



n = l,...,Ar-l, 
n = 2,...,N. 



(4.3a) 
(4.3b) 



By use of z^{n) and z^(n), we define the stationary-action path for the center position and 
the chirality as 



q\n):={z\n) + -z\n))/2, 
0«(n) := {z\n)-z''{n))/2i, 



(4.4a) 
(4.4b) 



which are not real in general. 

The left-hand side of ( [4.3| ) are 0{t) because of the factor e on the right-hand side. Hence, 
as far as the equations for the stationary- action path is concerned, we can go over to the 
continuous time: 



-insiTi\i{z\t)-z\t)), 
-zfisinh(^'^(t)-z^(t)). 



dz\t) 

dt 
dz%t) 

dt 



Note that the boundary condition is dictated by ( ^^ ) as 

z\0) = zi, r(T) = 4. 

It follows from Eqs. ( |4.5|) that 

d 



dt 



{-z\t)-z^{t)) = 0. 



Hence, 



z\t)-z\t) 



-2ic 



where (^ is a complex constant (0 = 0' + z0"; 0', (p" G R). Substituting (^78| 
taking account of the boundary condition (|4.6| ), we obtain 

^'(t) = -fitsin20 + zi, r{t) = -Q{t-T)sm2(l) + Zp. 



(4.5a) 
(4.5b) 

(4.6) 

(4.7) 

(4.8) 
into ( [4.5|) , and 

(4.9) 
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Putting this back into ( |4.^ ), we find that the constant is determined by the following 
algebraic equation: 



or equivalently, 



r2Tsin 2(1) + Zp — zi = — 2i0, 



20" — QT sin 20' cosh 20" = g = gp — gj, 
20' + nr cos 20' sinh 20" = 0f + 0i. 



(4.10) 



(4.11a) 
(4.11b) 



The stationary- act ion path can be expressed in terms of q^{t) and 0'^(t) as defined by 
as 



g^(t) = -Qt sin 20 + zi - 20, 



(4.12a) 
(4.12b) 



Note that both of these are complex. Eq. ( [4.12b| ) shows that (ff{n) — (jf{n — 1) = for all n 
. This justifies our procedure of working with ( |3.11D instead of (PTP|). Incidentally it follows 
from (|]T|) that 



dq^t) 
dt 



-fisin20"(t). 



(4.13) 



Thus, the velocity of the center position depends on the chirality and is proportional to the 
transverse anisotropy (recall that Vt (x. a). In the special circumstance that 0^(t) happens 
to be close to ±7r/2, we can put (jf{t) = ±7r/2 + (p^{t) to find 



dq%t) 
dt 



oc ip'it). 



(4.14) 



This shows that in such a circumstance the center position and the chirality are mutually 
canonically conjugate, thereby confirming the claim made in Ref. 10. In any case it is clear 
that q^{t) and 0'^(t) are closely coupled; they should be treated on an equal footing. 

B. Stationary action 



The stationary action iS'^fz'^, z^] may be arranged as 



ssc I cssd 



^ss .-s'[z',z'] =S'''' + S' 



n 



5^^^ := iVnw5 



l{\z^\' + \zr) + l{4AN-l) + z%l)zd 



N-l 



+ - E {(^'(^ + 1) - z\n))z\n) - z\n){z\n) - z'{n - 1))} 

2 n=l 






N 



"- n=l 



a 



1 + - {1 + cosh(r(n) - z\n - 1))} 



(4.15a) 

(4.15b) 
(4.15c) 
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The first and the second terms on the right-hand side of ( |4.15b|) depend on the initial 
and the final state. The second term depends also on T through the stationary-action 
path. These terms, which have been neglected in the continuous-time formalism, turn 
out to be crucial for a correct evaluation of the transition amplitude. The third term 
corresponds to the "Berry-phase term" in the continuous-time formalism. In the latter 
formalism, interference effects have been ascribed to this term alone. However, the second 
term can also contribute to interference effects. This is another remarkable difference from 
the continuous-time formalism. Of course, interference effects in question can arise only if 
there exist two or more stationary-action paths. In the case of a free domain wall under 
consideration, there is no question of interference because the stationary-action path is 
unique. 

Substituting ( ^^91 ) into ( [4.15| ), we obtain 



^ CSC 



NbwS 



I^fT + kiT - 24^1 - (firsin20)2 



Oie) 



-s 



sd 



h' 



EnwT 1 



- iNYiwSnr COS 20 + C(e) 



(4.16a) 
(4.16b) 



Putting these together, we finally arrive at 



h 



S' 



A^owS* 



\zf\^ + \zi\^ - 2z^zi - (fiTsin20)^ + 2ifiTcos20 



--£owr(i + - 



(4.17) 



This stationary action is complex in general. 



C. Fluctuations 



As noted previously,li3 well-defined evaluation of the fluctuation integral is possible only 
in the discrete-time formalism. Accordingly we separate the integration variables as 



z{n) =z''{n) + C{n), z* {n) = z"" (n) + C{n) : n = 1,2, ...,N - 1. 

It is convenient to define 

C(0) = C{N) = 0. 



(4.18a) 



(4.18b) 



Substituting ([4.18a|) into the action (|3.8dD and (|3.11| ), and expanding up to the second order 
in the fluctuation, we get 



S^[z\z]=S^^ + S2[CX], 



(4.19a) 
(4.19b) 



where 
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h 



S2[C, C] := -^^ E {2C(n)CH - 2C{n)an - 1)} 



(4.20a) 



n=l 



^ -dr/-* /•! .^OwS" 



^5,^[r,c]:= 



TV 

ficos20 E e {(CH)' + iCin - l)f - 2Cin)an - 1)} . (4.20b) 



n=l 



Accordingly, the transition amplitude reduces to 



{zp\e ''^'^^^\zi) ~ exp 



h 



S' 



K2{T), 



1 da n)dC{n) ft 



(4.21a) 
(4.21b) 



n=l 

where A^ is a constant, whose value as well as the detailed evaluation of K2{T) are given in 
Appendix B. The result is 



K2iT) 



Q^ = Q cos ; 



'I - in^T 
This completes a microscopic evaluation of transition amplitudes for a free domain wall. 



(4.22) 



V. TRANSITION PROBABILITY 

We can now compute transition probabilities between various initial and final domain- 
wall states: 



P(g,0F,0i;T):=|(^p|e-^^^/^|;^i)r 



gOn^T 



1 + 2'^Qd^T + IQJ^T^ 



exp 



-1 535^^ 
n 



h 



535^^ = A^Dw^ 



q" + {<f>F 



- 2(nT)2gfJ(sin20)2 - 2fiT3cos20 



(5.1a) 



(5.1b) 



The right-hand side does not depend on gp and qi separately but only on q{= qp — qi) as 
expected from the translation invariance. This justifies the notation P{q, 0f, 0i; T). 

A. Case of a = 

In the absence of transverse anisotropy (f? oc a = 0), the transition probability is inde- 
pendent of T: 



P(g, 0F, 0i; T) = exp [-NnwS {g' + (0f - 0i)'}] = | (^rki) T- 



(5.2) 



Thus, the transition probability coincides with the overlap between the initial and the final 
states and is a function of the differences gF — qi and 0f — 0i- It is shown in Fig. 3. 

The factor AdwS* in the exponent of ( p.2| ) is large (typically of order of 10^), reflecting 
the semi-macroscopic character of a typical domain wall of interest. Hence, if the final state 
is even slightly different from the initial state, the transition is forbidden. 
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B. Case of a / 

In the presence of the transverse anisotropy, the transition probabihty depends on T: 

P(g,0F,0i;T) 

exp(— r sin 20' sinh 20") 

1 - 2r sin 20' sinh 20" + r2{(cos 20' cosh 20") 2 + (sin 20' sinh 20") 2} 



X exp 



- A^DW>^{g' + (0F - 0l) 



-r^ [(sin 20' cosh 20")^ - (cos 20' sinh 2(f)")A + 2r sin 20' sinh 20"| 



(5.3) 



where r = QT. 



1. analytical evaluation in linear approximation 



In order to find 0(= 0' + i0"), we need to solve the algebraic equation ( [4.11| ). On 
inspection we see that it has a solution = (0f + 0i)/2 if g = — rsin(0F + 0i). This 
motivates us to look for a more general class of solutions by linearizing ( [4.11| ) under the 
following condition to be justified a posteriori: 



iF + 0I 



I ' J," " I 'I I "I ^^ 1 



(5.4) 



Then, we can write down the linearized version of ( [4.11 



as 



where 



Hence 



2(/?" - tCyi2^' = q', 2ip' + tCfi2^" = 0, 



g' = g + r^Fi, Cfi = cos(0F + 0i), ^fi = sin(0F + 0l) 



2lp' 



tQ 



FI , 



1 + (rCFi)2 



2ifi" 



1 + (rCFi)2 



q'. 



(5.5) 



(5.6) 



(5.7) 



Since the factors multiplying q' are at most of order unity, the assumed condition ( |5.4| ) is 
satisfied if \q'\ ^ 1. Substituting (|5.7| ) into (|5.1bD we find 



n 



J2 



+ 



.1 + (rCFi)2 
while the prefactor of (|5.1a| ) is obtained, up to 0{q''^), as 



(5.8) 



1 + (tCfi)^ 



tCfi 



1 + (tCfi)\ 



2TS-piq' 



2 \Cyi 



l + irC^iYT 



,'2 



(5.9) 
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Because of the large factor A^dw-S* in the exponent of ( p.8|) , the g'-dependence of the prefactor 
is neghgible. Hence, 



P{q,4>F,4>hT) 



1 



l + (rcos(0F + 0i))^ 
X exp 



^[ (g + rsin(0F + 0i))^ , 




(5.10) 



which is vahd provided that \q + r sin(0F + 4>i)\ ^ 1- This approximate formula shows that 
transitions are suppressed if 0f 7^ 0i- It also suggests the following picture: for a given 
(0F,0i), the peak of the wave packet representing the center position q moves with the 
velocity (in units of Af]) 

I'packet = - sin(0F + 0l), (5.11) 

while the width (in units of A) of the wave packet increases as 

Wpackct(r) = {1 + (rcos(0F + 0i))'}'/'. (5.12) 

The minimal velocity and maximal spreading occurs at |0f + 0i| = (mod it), while the 
maximal velocity and minimal spreading occurs at |(/)f + 0i| = tt/2 (mod tt). We choose 



these cases as well as the intermediate case of |0f + 0i| = 27r/3, and depict (|5.10| ) in Figs. 
4—6 with solid curves. These figures show the case of minimal suppression {(pp = (pi). We 
take NdyjS = 100 throughout (and also in Fig. 7 to be mentioned below). 

Fig. 4 shows the case of (pp = (pi = n/2, where fpacket = and Wpacket('7") = (1 + t'^Y^'^. 
The curve for q = 0, where the formula ( |5.10| ) is exact as mentioned at the beginning of 
this subsection, shows that the probability of remaining in the initial state decreases with 
time. Thus the domain wall exhibits a quantum phenomenon analogous to the wave-packet 
spreading for a free particle. The long-time tail originates from the prefactor of (|5.3|) coming 
from the fluctuation. If the final state corresponds to the mere displacement of the initial 
domain wall by the distance q ^ 0, the transition probability exhibits an initial increase, 
which may be interpreted as the appearance of an overlap between the final wave packet 
and the wave packet evolving from the initial state. For a small q, the overlap initially 
suppressed by the large factor A^dwS* is rapidly recovered resulting in a sharp initial increase 
of the transition probability. The origin of the long-time tail is the same in the case of 
q = 0. The case of (pp = (pi = -n/S, where f packet = ^3/2 and Wpacket('r) = (1 + rV4)^/^ is 
shown in Fig. 5, which exhibits a moving and spreading wave packet. So far, the quantum 
dynamics of the domain wall resembles that of a free particle. However, Fig. 6, where 
(pp = (pi = — vr/4, reveals quite a different feature; it may be interpreted as showing a non- 
spreading wave packet hke a solitary wave with f packet = 1 and Wpacket(T) = 1. Finally, Fig. 
7 shows that the transition probabilities between states with the common center position 
(g = 0) and different chiralities are rather small. This is in contrast to the situation in 
Fig. 4, where there is an appreciable transition probability between states with the common 
chirality {(pp = (pi) and g 7^ 0. 
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2. numerical evaluation 



To check the accuracy of the approximate formula ( |5.10| ), we have numerically solved 
( [4.11| ) for cf)' and 0", and put them into ( p.3|) . The results are depicted in Figs. 4—7. It is 
confirmed that ( [5. 101) remains valid even for \q + rsin(0F + 0i)| ~ 1- 



C. Comparison ^vith quantum dynamics of a free particle 



Let us compare the above-found quantum behavior of the domain wall with that of a 
free particle. 

The transition amplitude for a free particle with mass m between the initial coherent 
state \Zi)[Zi = Qi + iPi) and the final coherent state \Zp){Zf = Qp + iPp) is given by 



(ZF|e-*^^/"|Zi) 



lim 

i ^ 

tS[Z\Z]:=Y. 

'^ n=l I 



^^ dZ{n)dZ*{n) fz 

n=l 



27Ti 



h 



(5.13a) 



^{\Z{n)\' 



\Z(n-l)\' + Z*(n)Z(n-l] 



-zUl-{Z*{n)-Z{n-l)r} 



(5.13b) 



This action resembles that of a free domain wall ( |3.8d|) and (|3.11|) . Indeed, if we expand the 
non-linear part in ( p.8d| ) as 



cosh{z*{n) — z{n — 1)) ~ 1 



(5.14) 



then the action (|3.8| ) for a free domain wall reduces to ( ^.13b|) . In this sense, a free domain 
wall may be said to be a non-linear version of a free particle. 

The transition probability for the free particle is exactly calculated as 



X exp 




^+{i^y 



-(P, + fl)} -5ip(^.-fl: 



(5.15) 



where Qi and Qp are the center of the initial and the final Gaussian wave packet, respectively, 
and likewise Pi and Pp are the initial and the final mean momentum. The width of the wave 
packet has been chosen to be 6 both for the initial and the final state. The form of (|5.10|) 
differs from that of ( |5.15| ) only in the non-linear factors cos(0f + (pi) and sin(0F + 4>i)- If 
these are linearized as cos(0f + 0i) ~ 1 and sin(0F + (pi) — (pF + 4>i^ then (|5.1CI|) reduces 
to the same form as ( 5.15 ); in this case the behavior of a free domain wall is the same as 
that of a free particle. However, the linearization is not always allowed. In particular, when 
(pF + <Pi = i'?'"/2, the behavior of the domain wall is completely different from that of the 
free particle; the wave packet does not spread! This is a manifestation of the non-linear 
character of spin. 
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D. Effective mass of a free domain wall 

We can estimate the effective domain-wall mass from the correspondence of the domain 
wall and the particle as noted in the previous subsection. 

Comparison of ( |5.10| ) and ( |5.15| ) reveals the following correspondence. First, the coeffi- 



cient of T in the prefactor suggests 



fi|cos(0F + 0i)|^-^^. (5.16) 



Second, the coefficient of q^ in the exponent at T = suggests 

A 



JDW 



, = < — > 6. (5.17) 

(5dw can thus be interpreted as the width of the wave packet describing the initial domain 
wall. These two correspondences suggest to associate the domain wall with the "effective 
mass" Mdw given by 

°^ 2fi|cOs(0F + 0l)|^DW |cOs(0F + 0l)r 

It is to be noted, however, that this "effective mass" depends on the initial and the final 
chirality; as such, it can not be viewed as an effective mass of an ordinary dynamical entity. 
It coincides with the Doring mass if 0f + 0i = or ivr. On the other hand, it is infinite 
if (pF + 4>i = ±7r/2, which is just another way of expressing the non-spreading of the wave 
packet as noted in the previous subsection. 

VI. DISCUSSION 



If we calculated the transition probability by use of the continuous-time action ( |3.12a| ) 



what result would have been obtained? Though the stationary-action path is the same as 
^, the stationary action would be 



- -EdwT (^1 + - j . (6.1) 

(6.2) 

'cos20'sinh20" - sin20'sinh20"). (6.3) 

Let T = 0. Assuming that the fluctuation integral K2con{T) was somehow evaluated and 
that it agreed with the correct value ( [4.22 ), we would have found \{z-p\e~'^^'^^^\ziiW^Q ~ l- 





^ CSS . 


- '9^ 
^'^con 


[r,z'] 






= 


= —iND^jv St {(f) sin 2(f) 


+ cos 


Hence 










1 (^F e" 


^HT/n^^^^2 


~ |K2con(T)2exp('- 




-Us 

h 


'con = ^DW 


.5r(0"sin20'cosh20' 


' + <P' 



However, this is a completely meaningless result, since the left-hand side should be equal 
to I (zf 1 2:1) p. (If one started from the effective action ( |3.15| ), one would also flnd a similarly 
meaningless result.) Furthermore one can not rationally calculate the fluctuation integral 
from this formalism.lill 
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VII. CONCLUDING REMARKS 

We have considered the macroscopic quantum dynamics of a free domain wall in a 
quasi-one-dimensional ferromagnet by use of the spin-coherent-state path integral in the 
discrete-time formalism. The center position and the chirality, which have been chosen as 
the collective degrees of freedom, are noted to be mutually canonically conjugate in a loose 
sense. The quantum behavior of the domain wall is the same as that of the free particle and 
its effective mass is the Doring mass if 0f + 0i = or ±7r, but in general it differs from the 
latter in some non-linear effects. We have also pointed out some grave difficulties associated 
with the continuous-time formalism. It can not correctly evaluate transition amplitudes. 
Its assertion on interference effects on the basis of the "Berry-phase term" alone is also 
questionable. 

Let us speculate on MQT and MQC involving a domain wall. Since a free domain 
wall with a fixed chirality ((/>f = 4>i) has been shown to behave roughly like a free particle 
unless (p-p = (pi = ±7r/4, we expect that a quantum depinning (MQT) will occur in the 
case of a weak pinning and a strong transverse anisotropy as mentioned by many workers; 
a strong transverse anisotropy tends to fix the chirality at 0f = 0i = ±vr/2. However, if 
a transverse anisotropy energy is comparable to a pinning potential, the dependence of the 
domain-wall mass on the chirality can be important. This may somewhat affect the MQT. 
Such a possibility has been overlooked in the literature. The MQC has been suggested to 
occur in the case of a strong pinning and a weak transverse anisotropy,^ namely for a fixed 
center position (g = 0). However, Fig. 7 shows that the transition probability for q = and 
I0F — (^i| ~ TT is negligible. This is due to the large factor A^dwS* in the exponent. Hence, 
for the MQC to occur, it may be necessary to invoke a mechanism (e.g., a magnetic fieldB) 
to decrease {(pp — (f)i\. At any rate, a careful consideration is needed to make a conclusion 
on the possibility of the MQC. 
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APPENDIX A: DERIVATION OF EQ. (3.9) 

The numerator in the logarithm of ( p.8c|) can be rewritten in two forms; 

-^ _|_ ^-x+z*{n)+z(n-l) 

I p-x-\-z*{n)^'z{ri) \ 

= (1 + e-+-*W+-W) 1 + — ^, , , , (^-{An)-zin-X)) _ ^) ,p^^. 

^ ' \ \ j^ g-x+2*(n)+2;(n) ^ ' y ^ ' 




]^ _j_ p—x^z*{n—\)^z{n—\) 

In order to perform the x integration, we pay attention to the following formula; 



da; In I 1 



A 



1 + e^'+^ 



Di 1 



A 



Di 



■1 + 



l+e^+-B 



+ ln(l + A)ln 



1 + e^+^y "V 1 + ^ 



1 + e^+^ / ' 



where Di(z) is the dilogarithm defined bycj 



Diiz) :-- 



dt 



Int 



'i 1-t 
Thus, the integral in ( p.8c|) can be evaluated as 



/ 



L/A 



(ixln 



(1 _|_ ^-x+z*{n)+z{n-l)\2 



j^j-y l\ _j_ g-a:+z*(n)+^(n)V]_ _j_ ^-x+z* {n-\)+z(n-\)\ 



E 

Q = l 



/^? 1 + 



Ar 



\ _|_ gX+Ba 



Di 



1 + r 



i + A 



+ ln(l + A^) In 






L/A 



J -L/A 



where 



J^ = g-{2W-^("-l)) _ ]_ = g-('7(")-'?{"-l))-J('/'o(n)-</'o(n-l)) _ ]_ 

5i = -{z\n) - z{n)) = -2q{n), 

A ^ g2*(n)-2*(n-l) _ j^ _ g9(7i)-g(n-l)-i(0o(n)-(/io(n-l)) _ ]_ 

^2 = -(2*(n - 1) - z{n - 1)) = -2q{n - 1). 
Since L/A ^ 1, (|A4|) may be simplified as 



a=l 



Di{l + A^) + Dt ( Y^) } + 1^(1 + ^") (x + l^exp(-5„ 



Next, we use the dilogarithm identity 

Di{A) + Di 



A 



-(InAf 



to find 

2 



o=l L 



(A2) 



(A3) 



(A4) 



(A5a) 
(A5b) 
(A5c) 
(A5d) 



(A6) 



(A7) 



where 



^(ln(l + A^)f + ln(l + A^) (^j + Inexp(-Bj) 
[(gH - qin - l))^ + i?(A0oH) + t {2(gH + g(n - 1)) + 2L/A} /(A0oH) 



A0o(?^) = M^) - M^ - 1). 

R{(f)) ■= - (lnexp(-#))^ : |0| < vr, 

/(0) := ilnexp(— i0) : —n < (p < n. 



(A8) 



(A9a) 
(A9b) 
(A9c) 
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R{4>) and /(0), which are 27r-periodic functions, can be cast into the form (|3.10| ). This 
periodicity follows from (|3.8c| ), which originates from the overlap of between spin-coherent 
states ( p.9| ) whose real and imaginary parts are even and odd periodic, respectively. The 
discontinuity or non-smoothness are the consequence of the spatial continuum approxima- 
tion. 



APPENDIX B: CALCULATION OF FLUCTUATION INTEGRAL 



The fluctuation action ( |4.2(]| ) may be written as 



-S-{N-l):=-z- 



h 



Y.{MnW{n)f + B{n){an)f 



n=l 



where 



+2C(n)C(n)CH + 2D{n)Cin)an - 1)} , 



A{n) = £^2 cos 20, 
B{n) = eQcos2(f), 
C{n) = -I, 
D{n) = i — eQ cos 2(f). 



This can be expressed in the matrix form: 



where 



^SnN-l) = -^^'CM{N-l)C 



^C := (r(iV - 1), C(iV - 1), C(iV - 2), ..., C(l), C(l)) 



(Bl) 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



(B3) 



(B4a) 



M{N-1) :-- 



(B4b) 



fA{N~l) C{N-1) \ 

C{N-1) B{N-1) D{N-1) ■•• 

D{N-l) A{N~2) ■•. 

••. D{2) 
D{2) A{1) 
V C{l)B{l)) 

As to the integration measure, which has not been explicitly derived, we may make the 
following two assumptions: 

(i) Its structure is the same as that of the (boson-) coherent-state path integral: 



1 dan)dC{n) 



M 



27X1 



(B5) 



where M. is a. constant. 
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(ii) The value of Ai can be inferred from the overlap between the initial and the final 
states. 

Assumption (ii) is reasonable since exp[iiS'^'^(T = 0)/h] should coincide with {zp\zi). On 
the basis of these assumptions, the fluctuation integral can be cast into the form 






1 dC{n)dC{n) 



nil M 



2m 



exp 

\N~-1 



-z^^^^*CM(iV-l)C 



lim UMNj,^Sf^''-'\-l)^~-' det{iM{N - 1)) 



TV^oo 



-1/2 



lim 



{MN^wSf^^-^'^ det M{N - 1) 
The above determinant can be evaluates as follows.ll3 Let 



-1/2 



(B6) 



M(n) := det M(n) 



A{n) -i 

-i B{n) D{n) 
D{n) A{n - 1) 



(B7) 



This can be expanded in terms of the cofactor as 

M{n) = A{n)M'{n) + M(n - 1), 
where 



(B8) 



M'(n) := 



Bin) Din) 

Din) Ain- 1) -i 
-i Bin -I) 



(B9) 



M'in) can in turn be expanded as 

M'in) = Bin)Min - 1) - D^in)M'in - 1). (BIO) 

The recursion relations (|B8D and (|B1CI|) should be solved with the initial condition 

M(0) = 1, M'(0) = 0. (Bll) 

In the limit of e — > 0, they reduce to a set of coupled first-order differential equations: 

dMit) 



dt 
dM'it) 



%M'it), 
--n^{Mit) + 2iM'it)}, 



(B12a) 
(B12b) 



M(0) = 1, M'(0) = 0, 
where Q^ is given by ( [4.221) . This gives 
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ii L(kt 



Hence, 



M{t) = (1 - in^t)e'^^^\ M'it) = n^te'^^^\ 



K2iT)=\im {MNBwSf^^-^^M{N-l] 



-1/2 



g-iQ^T/2 



lim 



Since -^'2(0) should be unity, we conclude from assumption (ii) that 

1 



M 



NnwS- 



(B13) 



(B14) 



(B15) 
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FIGURES 

FIG. 1. Domain walls with three chiralities (quoted from Ref. 10); (a) right-handed wall 
((/>0 = 7r/2), (b)left-handed wall ((/)q = —tt/2), and (c) wall with no chirality (0o = 0). Circles in 
(a) and (b) drawn to guide the eye lie in the yz plane, while the spins lie in the zx plane in (c). 
The quasi-one-dimensional direction of the crystal is here aligned with the spin hard axis for ease 
of visualization. A different alignment, which may be the case for a real magnet, does not affect 
the content of the text; for instance, one could rotate all the spins by tt/2 around the y axis if the 
dominant anisotropy originates from the demagnetizing field. 



FIG. 2. Functions R{4>) and /(0). 

FIG. 3. The transition probability in the absence of transverse anisotropy as a function of 
Q ^ Qf — Qi and 0f — 4'i with iVow'S' = 100. It is independent of T. 

FIG. 4. Time-dependence of the transition probabilities for (pp = cpi = tx 12. Solid lines represent 
analytical results in the linear approximation (the curve for g = is exact). Squares, open circles, 
and triangles represent numerical results. 

FIG. 5. The transition probability as a function of q and T in the case of 0f = 0i = — tt/S. 
Numerical results are indistinguishable from the analytical ones. 

FIG. 6. The transition probability as a function of q and T in the case of 0f = 01 = — '7r/4. 
Numerical results are indistinguishable from the analytical ones. 

FIG. 7. Time-dependence of the transition probabilities for gp = Q'l and ^i = 7r/2. Solid lines 
represent analytical results in the linear approximation. Squares represent numerical results. 
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